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We propose a simple theoretical model for a molecular chemical engine that catalyzes a 
chemical reaction and converts the free energy released by the reaction into mechanical work. 
Binding and unbinding processes of reactant and product molecules to and from the engine are 
explicitly taken into account. The work delivered by the engine is calculated analytically for 
infinitely slow ("pseudo-static") processes, which can be reversible (quasi-static) or irreversible, 
controlled by an external agent. It is shown that the work larger than the maximum value 
limited by the second law of thermodynamics can be obtained in a single cycle of operation by 
chance, although the statistical average of the work never exceeds this limit and the maximum 
work is delivered if the process is reversible. The mechanism of the energy transduction is also 
discussed. 

KEYWORDS: molecular chemical engine, quasi-static process, the second law of thermodynamics, energy 
transduction 



1. Introduction 

Biological molecular motors, such as myosin, kinesin, 
dynein and Fi-ATPase, are tiny machines that catalyze 
the hydrolysis reaction of adenosine triphosphate (ATP) 
and convert the free energy released by the reaction into 
mechanical work.^"^ The mechanisms of energy transduc- 
tion in these molecular chemical engines are not fully un- 
derstood. A simple, idealized model of a molecular chem- 
ical engine, like the Carnot cycle of a heat engine, would 
provide clear understanding of the energy transduction 
in biomolecular motors and insights into the design of ar- 
tificial nanomachines. The purpose of the present paper 
is to propose such a model. 

The theoretical models that have been proposed to 
study how the chemical free energy can be converted into 
mechanical work may be grouped into two classes. One 
class of models'''^ introduce a "reaction (or chemical) 
coordinate" x as well as a "mechanical coordinate" y. 
The former coordinate describes how the chemical reac- 
tion proceeds (the increase of x by, say, a will correspond 
to the consumption of one fuel molecule due to the reac- 
tion) and the latter indicates how the engine translocates 
along a linear track (with a periodic structure of period 
b) or how the rotor in the engine rotates (6 = 2tt). The 
"mechano-chemical coupling" is expressed as a potential 
function U{x,y) that satisfies U{x-\-a,y) — U{x,y)-\- AG 
and U{x,y -\- b) = U{x,y), where AG < is the change 
in the free energy of the reservoir containing fuel and 
product molecules due to the reaction of a single fuel 
molecule. In this approach, the binding of a fuel molecule 
to and the unbinding of a product molecule from the en- 
gine are not taken account explicitly. As a result, the 
reservoir comes into play only through AG, which is 
unsatisfactory because the chemical potentials (or the 
densities) of fuel and product molecules can affect inde- 
pendently the performance of the engine, as one might 
expect or we shall see in this paper. 



The other class of models deal with transitions be- 
tween discrete chemical states j {j = 1,2, ... ,n). This 
class can be divided into two subclasses according to 
whether the mechanical coordinate is discretized^ or con- 
tinuous.^"^" In the former, the value of the mechanical 
coordinate is determined uniquely by the chemical state, 
and hence a transition results in a translation of the en- 
gine (motor) or a rotation of the rotor. In the latter, 
the engine in state j undergoes Brownian motion under 
the influence of state-dependent potential Uj{y) associ- 
ated with the mechanical coordinate y, and the transi- 
tion from one state to another does no accompany the 
translational or rotational motion of the engine at the 
instance of transition. In this class of models, the bind- 
ing of a reactant molecule to the engine, the chemical 
reaction, or the unbinding of a product molecule from 
the engine trigger a transition between particular states. 
However, the mechanism of energy transduction is not 
quite transparent in the previous works on these models. 

In the present work we take a combined approach: a 
continuous reaction coordinate is used to describe the 
chemical reaction catalyzed by the engine, and discrete 
chemical states are introduced to distinguish the state of 
the engine with a fuel or product molecule bound to it 
from the state with no bound molecule, while a continu- 
ous mechanical coordinate denotes the rotation of a rotor 
connected to the engine. In constructing the model pro- 
posed here, we are inspired by the models for Fi-ATPase, 
a biological rotary motor, by Wang and Oster^ and by 
Kinoshita et al.^" This protein machine is remarkable 
in that it works not only as a motor^-'^'^^ (when ATP 
molecules are consumed) but also as a generator^^' 
(ATP molecules are generated when it is forced to ro- 
tate in the reverse direction), just as a heat engine works 
as a heat pump when it is run in the reverse direction. 
The minimal ingredients necessary to build a reversible 
molecular chemical engine are extracted from their mod- 
els®' and an assumption on the chemical reaction is 
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introduced to simplify the analysis. 

Our model is simple enough to analytically work out 
extremely slow processes controlled externally, which are 
analogous to the Carnot cycle of a heat engine. We will 
call such a process a pseudo-static process instead of a 
quasi-static process, because it turns out that such a slow 
process is not always reversible in the thermodynamic 
sense. It will be shown that the average work delivered 
by the engine is identical to the free energy (—AG > 0) 
released by the chemical reaction if the process is re- 
versible, and it is less than |AG| if the process is irre- 
versible; this result is consistent with the laws of ther- 
modynamics. However, it is found that in a single cycle 
of operation, work larger than the maximum work im- 
posed by thermodynamics can be obtained by chance if 
the process is irreversible. 

In the present model it is the torque resulting from 
the interaction between the engine and a bound (fuel or 
product) molecule that contributes to the net mechan- 
ical work delivered in a cycle of pseudo-static process, 
and this torque is independent of the chemical potentials 
of fuel and product molecules in the reservoir (i.e., inde- 
pendent of AG). It will become clear that the chemical 
potentials determine the "timings" of the binding and 
unbinding of fuel and product molecules, thereby con- 
trolling the interval over which the torque is in effect. 
Thus the amount of work delivered by the engine is af- 
fected by the densities of fuel and product molecules in 
the reservoir. 

2. Model 

We wish to construct a simple model for a molecular 
chemical engine that is supposed to utilize a hypothetical 
chemical reaction 

A^B. (1) 

The setting we assume for the engine to work is as fol- 
lows. The engine is coupled to a rotor, whose rotational 
angle 9 can be controlled by an external agent. The en- 
gine exchanges molecules A and B with a reservoir that 
contains a large number of these molecules. The chemical 
potentials ha and /is of molecules A and B, respectively, 
are assumed to satisfy the condition 

All = fj,A - IJ-B > 0. (2) 

The engine and the reservoir is in contact with a heat 
bath of temperature T. The chemical reaction (1) under- 
goes only at a "binding site" on the engine, to which a 
single molecule A or B can be bound, and the engine has 
only one binding site. 

It will be said that the engine is operated forward or 
backward depending on whether the angle 6 of the rotor 
is increasing or decreasing. The engine is designed so that 
it cannot bind a molecule when ^ = (mod 27r). A for- 
ward (backward) cycle is a process in which 9 is increased 
(decreased) from zero to 27r (— 27r); the engine is empty 
both at the start and the end of a cycle. A forward cycle 
is expected to proceed as follows. A "fuel" molecule A is 
bound to the engine. Then a "forward reaction" A ^ B 
takes place at the binding site. Finally, a "waste" (prod- 
uct) molecule B is released to the reservoir. Meanwhile, 
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the engine delivers work by exerting torque on the rotor. 
If the cycle proceeds as expected, the free energy of the 
reservoir changes by AG = —Afj, < and the engine 
returns to the original state after the cycle. If the laws 
of thermodynamics can be applied to this process, the 
maximum work the engine can deliver in principle will 
be A/X.15 

Now we describe in detail how our model is con- 
structed. 

2. 1 Liga,nds 

Molecules A and B are thought of as two stable states 
of a single molecule, which will be referred to as a ligand, 
with different geometrical structures. It is assumed that a 
continuous transformation of structure from A to B can 
be characterized by a single, dimensionless parameter r, 
a "reaction coordinate." 

The energy profile associated with the structural 
change of a ligand molecule U{r) may be defined, apart 
from an unimportant additive constant, by 

exp [-/3[/(r)] = Tr cxp [-^Hl (r)] , (3) 

where f3 = l/ksT with fee being the Boltzmann constant 
and Hl (r) is the Hamiltonian for the vibrational and ro- 
tational modes of a ligand molecule under the constraint 
that the reaction coordinate r is held fixed; we assume 
that r can be treated as a classical dynamical variable. In 
the definition (3) of U{r), we have assumed, for simplic- 
ity, that the reservoir is an ideal gas mixture of molecules 
A and B. If the reservoir is a dilute solution where A 
and B are solutes, the interaction of a ligand with sol- 
vent molecules should be taken into account in defining 
U{r). Once U{r) is defined appropriately, the discussions 
presented below will be valid irrespective of whether the 
reservoir is an ideal gas or a dilute solution. 

In order for molecules A and B to be stable, U{r) must 
have two minima, say aXr = rA and r = rB, correspond- 
ing to molecules A and B. The height of potential barrier 
separating these minima should be much larger than the 
thermal energy kuT so that reaction (1) does not occur 
spontaneously in the reservoir. 

Let riA and n_B be the densities (numbers of molecules 
per unit volume) of molecules A and B, respectively, in 
the reservoir. Then the chemical potentials of A and B 
are given by 

n 

/x„ = C/„ + fcBrin— {a = A,B), (4) 

where Ua = U{ra), and rja is a quantity that has di- 
mensions of density and is proportional to a restricted 
partition function associated with the reaction coordi- 
nate: 

r/„ oc / exp {-/3[f/(r) - f7„]} dr {a = A, B). (5) 

Here, TZa is a small interval around r = on the r axis; 

if U{r) can be approximated by U{r) = Ua + \Ka{r — 
rocY with a positive constant in the vicinity of ra, 
we will have rja oc (27r/3/ii'a)^/^. 
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Fig. 1. The conformational change of the engine upon the rota- 
tion of the rotor is schematically shown. The bold arrows indicate 
the direction of a forward operation, (a) 6 ~ 0, (b) < 9 < Si , 

(c) e ~ 01, (d) 01 < e < 02, (e) e ~ 02, (f) 02 < < 2tt. 

2.2 Engine with no ligand 

As already explained, the engine is coupled to a rotor 
(through an arm- like structure, for example). We imag- 
ine that the engine is a macromolecule like a biological 
molecular motor and assume that its structure (confor- 
mation) changes with the rotation of the rotor. In this 
subsection we consider the case that the engine carries 
no ligand molecule. 

Let Hq{9) be the Hamiltonian of the engine with no 
ligand under the condition that the rotation angle of the 
rotor remains fixed to be 9. The potential (free energy) 
Vo(^) of the engine under this condition is defined by 

exp [-mm = Tr exp [~(3Hoie)] . (6) 

It is assumed that the engine is operated {9 is varied by 
the external agent) so slowly that it is in thermal equi- 
librium at any instant. Then the torque to{9) exerted by 
the engine on the rotor is given as the statistical average 
of —dHo{9)/d9, which imphes 

. (7, 

Obviously, Vo{9) and to{9) possess the following period- 
icity: 

Vo{9 + 27:)^Vo{9), Toi9 + 2n) ^ to{9) . (8) 

In order that the engine works properly and can be 
analyzed without much difficulty we require the engine 
to possess the following properties associated with its 
conformational change. Let 9i and 92 be some constants 
satisfying < 6'i < 02 < 2tt. When the angle 9 of the 
rotor is in the interval < 9 < 61, only molecule A 
can be attached to or detached from the binding site of 
the engine; if < 6* < 6*2, no ligand can be attached 
or detached; and when 92 < 9 < 2tt, only molecule B 
can be attached or detached. This can be made possible 
if, for example, the engine has a window-like structure 
and a ligand can reach the binding site only by going 
through the opening of the window (see Fig. 1): when 
< 9 < 9i, the opening of the window has a shape such 



Fig. 2. Potential W{9,r) defined by eq. (9) as a function of r for 
fixed values of 0. (a) ~ 0, (b) < < 0i, (c) ~ 0i, (d) 
01 < < 02, (e) ~ 02, (f) 02 < < 27r. 

that a molecule A can go through it but a molecule B 
cannot; the window is closed if 6*1 < 6* < 6*2 or 9 0; 
and when 6*2 < 9 < 2tt, the shape of the opening allows 
a molecule B but not a molecule A to go through it. 

2.3 Engine with a ligand 

Now we consider the case that a ligand molecule is 
bound to the engine. Let Hi-at{9^r) be the interaction 
between the engine and the ligand under the constraint 
that the angle 9 of the rotor and the reaction coordinate 
r of the ligand are held fixed. We define potential W{9, r) 

by 

exp{-/3[K)(0)+W^(0,O]} 

= Tr exp{-/3[i/o(f?) + i?L(r) + ffi„t(0, r)]} , (9) 

where Vq{9) has been defined in eq. (6). Roughly speak- 
ing, W is the sum of the ligand energy and the engine- 
ligand interaction. 

The performance of the engine depends crucially on 
potential W{9, r). Let us describe the desired dependence 
of W{9,r) on r for fixed values of 9 (see Fig. 2). If the 
angle of the rotor is in interval {) < 9 < 9i^W should have 
a deep well around r = so that a molecule A sitting at 
the binding site stays there as molecule A without being 
converted to molecule B] it may have another minimum 
around r = rs separated by a large barrier. Similarly, we 
require that W should have a deep well around r — tb if 
92 < 9 < 27r, so that molecule B is stable at the binding 
site. A ligand bound to the engine is supposed to undergo 
forward reaction A ^ B if 9 is varied from 9i to 6*2 or 
backward reaction B ^ A if 9 is varied from 92 to 9i . In 
this regard it is demanded that W has only one minimum 
around r = rA for 9 9i, and only one minimum around 
r — rs for ^ 02. For 9i < 9 < 92, W may have a single 
minimum at some r between and rs , and its location 
moves continuously as 9 is varied between 9i and 92] or 
W may have a secondary minimum slightly above the 
primary one separated by a barrier whose hight is much 
less than A;bT. 

For the sake of simplicity the relaxation of the reaction 
coordinate r is assumed to occur quickly enough that 
this degree of freedom is always in thermal equilibrium 
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during an operation of the engine. Then we can think of 
an "effective interaction" e{9) defined by 



exp I 



(3e{e)]= [ cxp[-(3W{9,r)] dr, 
Jn{e) 



(10) 



where the range TZ of integral depends on 6 as follows: 
TZ is the entire range of r for < 9 < 62, TZ = TZa 
[small interval around r = ta introduced in eq. (5)] for 
< 61 < Si, and 7^ = TZb for 02 < < 2n. The reason 
why the integral range is restricted to TZa for < 6 < 61 
is that as a function of r may have a valley around 
r = rs as mentioned above but this valley cannot be 
visited due to the large barrier separating it from the 
valley around r — rA- For the similar reason the integral 
range is restricted to TZb for 62 < 6 < Stt. Because of 
these restrictions on the integral range in eq. (10), e{0) 
is discontinuous at = (mod 2tt) in general, thought 
it is periodic, e{9 + 27r) = e{9). It is not difhcult to see 
that the engine carrying a ligand exerts torque 



nie) 



(11) 



d0 de 

on the rotor when its rotational angle is 9; it is the sta- 
tistical average of —d[Ho{9) + Hint{9,r)]/d9. 

2.4 Binding and unbinding of a ligand 

Suppose that the rotational angle 9 of the rotor is in 
range < 9 < 9i. Then the binding site of the engine can 
be empty or occupied by a molecule A. If the binding- 
unbinding process is in equilibrium, the probability pa(^) 
that a molecule A is bound to the engine is given by 

Zi{9)exp{f3fiA) 



Pa{0) = 



(12) 



Zo{9) + Zi{9) exp(/3/i^)^ 

where and Zi are, respectively, the partition func- 
tions of the engine without and with a ligand. From the 
definitions of Vb(0) and e{9) we have 

Zo(^^) =exp[-/3yo(e)], (13) 
Zi(0)=exp{-/3[K,(0)+e(0)]}, (W) 
and therefore eq. (12) can be expressed as 
1 



^^^^^ exp{/3[e(0)-AiA]} + l 
Similarly, we have 

1 



{Q<9<9i). (15) 



h<e < 2tt) (16) 



^''^^^ exp{/3[e(0)-A*B]} + l 

for the probability that a molecule B is bound to the 
engine when 9 is in interval 92 < 9 < 2TT. 

We have assumed that the engine is designed so that 
a ligand cannot be bound when 9 — (mod 27r). Consid- 
ering eqs. (15) and (16), this requirement implies that 

e^p{-P[e{+0)-fiA]}^0, (17) 

exp{-/3[e(27r-0)-^B]} = 0. (18) 

(Remember that in general e{9) is discontinuous at 9 — 0, 
i.e., e(0+) 7^ e(27r — 0).) In reality these conditions cannot 
be satisfied exactly, because e cannot be infinite. We con- 
sider eqs. (17) and (18) as mathematical idealization for a 
situation in which the left-hand sides of these equations 
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Fig. 3. An example of effective interaction e{6) between the en- 
gine and a bound ligand defined by eq. (fO). Binding and unbind- 
ing transitions of a molecule A occurs for < 8 < 81. Chemical 
reaction A <=^ B undergoes if a ligand is bound for 61 < 8 < 82; 
in this interval both binding and unbinding of a ligand is pro- 
hibited. Binding and unbinding transitions of a molecule B take 
place for 82 < 9 <2it. 



are extremely small and can be treated as zero practi- 
cally. 

Figure 3 shows an example of "effective interaction" 
e{9) between the engine and a ligand that satisfies con- 
ditions (17) and (18) and certain properties we are go- 
ing to discuss. When the engine is operated forward, it 
is preferable that probability pa{0) increases as 9 is in- 
creased from zero to 0i, while pb{9) decreases as 9 is 
varied from 92 to 27r. Hence, we assume that e{9) de- 
creases with increasing 9 iov Q < 9 < 9i and increases 
with 9 for 92 < 9 < 2tt. The dependence of e on 9 for 
9i < 9 < 92 can be rather arbitrary, because binding or 
unbinding of a ligand is prohibited in this interval. 

If the chemical potential fiA of molecule A lies well 
above ei = e{9i), the engine carries a molecule A with 
large probability when 9 reaches 9i in a forward opera- 
tion. Similarly, the engine carries a molecule B with large 
probability when 9 is decreased to 92 in a backward op- 
eration if /ZB hes well above £2 = £(^'2). These conditions 
are favorable for good performance of the engine. We 
shall not assume, however, that these are always satis- 
fied, because we are interested in how the engine may 
work under poor conditions as well as ideal ones. 

3. Pseudo-Static Operations of the Engine 

Let us analyze the work delivered by the engine when 
it is operated slowly enough that the binding and unbind- 
ing of ligands is in equilibrium at almost any instance of 
the operation. We call this process a pseudo-static pro- 
cess to distinguish it from a quasi-static process; the lat- 
ter is preserved for a reversible process, while the for- 
mer may be reversible or irreversible as we shall see. In 
a pseudo-static process the probability that the engine 
carries a ligand is given by eq. (15) for Q < 9 < 9i and 
by eq. (16) for 92 < 9 < 2tt, and the torque exerted on 
the rotor by the engine is given by 



r{9)^[l-PA{9)]To{9)+pA{9)n{9) 



(19) 



for < 6* < 01 and by 

t{9) = [1 - pb{9)]to{9) +pb{0)ti{9) (20) 
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for 02 < < 27r, where tq and ti are defined in eqs. (7) 
and (11), respectively. These expressions for the torque 
can be rewritten as 



(21) 



where a = A for < 9 < di and a = B for 02 < 
9 < 2n. Note that the quantity in the braces in eq. (21) 
is identical to the grand potential (grand-canonical free 
energy) of the engine. For 6i < 9 < 92 the torque 
is given by t{0) = to{9) or t{0) = t\{9) depending on 
whether the engine carries no ligand or one. 

3.1 Forward operation 

In a forward, pseudo-static operation of the engine the 
external agent is supposed to control the angle of the 
rotor to increase extremely slowly. The work delivered 
by the engine in one cycle of forward operation is given 

by 

/■27r 

Wfw= / T{9)de, (22) 



JO 



where t{9) is the torque exerted by the engine. 

We divide one cycle of a forward operation into three 
phases: the period in which < 9 < 9i will be called 
"phase A," the period oi 9i < 9 < 92 "phase C," and the 
last one {62 < 9 < 27r) "phase B." In phase A the torque 
exerted by the engine is given by eq. (21) with a = A, 
and therefore the work delivered in this phase is found 
to be 



Wft = Vo{0) ~ Vo{9,) - iln(l -pi) 



where probability pi is defined by 
Pi =Pa{9i). 



(23) 



(24) 



Condition (17) has been used in deriving eq. (23). The 
work obtained in phase B can be worked out sim- 
ilarly. Using expression (21) with a = B and condition 
(18), we get 

= Vo{92) - K)(27r) + 4 ln(l - P2), (25) 



where probability p2 is defined by 
P2 =Pb{92)- 



(26) 



In phase C there are two alternatives: if the engine carries 
no ligand, the work delivered in this phase is 



Vq{9^)-V^(92) 



(27) 



since the torque is given by tq {9) : and if a ligand is bound 
to the engine, the work is given by 

rCl 



fw 



Vo{9i)-Vo{92) + ei9i)-ei92 



(28) 



since the torque is ti{9). Summing up the contributions 
from the three phases, we obtain 



Wg 



1, 1-P2 



(29) 



for the work delivered in a cycle in the case that the 
engine carries no ligand in phase C. If the engine carries 
a ligand in phase C, the work in one cycle is given by 



(30) 



Let us examine the results for W^^ and Wj^ obtained 
above. It is noted that depending on the design of the 
engine and/or the preparation of the reservoir any of the 
cases (i) pi = P2, (ii) Pi < P2, and (iii) pi > P2 can be 
realized. In the first case we have 



(31) 



This consequence is in accordance with our naive expec- 
tation that the engine delivers no work if it does not 
catalyze the chemical reaction, and the free energy A/i 
released in reaction A ^ B is converted into mechanical 
work if the engine catalyzes the reaction. In the case of 
Pi < P2 one finds 



Wl<0, Wg>A^^. 



(32) 



Note that the second inequality tells that work greater 
than the the maximum work dictated by the second law 
of thermodynamics can be extracted in one cycle of op- 
eration. Apparent violation of the second law may occur 
also in the case of pi > P2, in which we observe 



> 0, 



wl < Am. 



(33) 



The first inequality implies that the heat extracted from 
a single heat bath can be converted into mechanical work 
without leaving any changes in the engine and the reser- 
voir of ligand molecules, which appears to be in contra- 
diction to the second law. 

The fact that the engine can deliver work larger than 
what is expected from thermodynamics in one cycle of 
operation does not necessarily mean inadequacy of the 
present model. This is because such an event occurs 
stochastically and in an uncontrollable manner. In the 
case of pi < P2, for example, the work delivered may 
be larger than the maximum work limited by thermody- 
namics {Wg > A/i) in one cycle but it may be less than 
the maximum work {Wg < 0) in another. Therefore, the 
appropriate quantity to be worked out should be the av- 
erage work per cycle in a forward operation carried out 
over many cycles. In a forward operation the probability 
that a ligand is bound in phase C is given by pi. Hence, 
the average work delivered per cycle is 

W^fw = {l-pi)W^^+piWg =piA/x- ■^/(pi,p2), (34) 

where function f{x,y) is defined by 



1 



f{x,y) = (1 -a;)lnY 



— -|- X In — . 

y V 



(35) 



It is easy to see that f{x,y) > for < x < 1 and 
< y < 1, where the equality holds ii x = y. Thus we 
find 



Wfw < Pi All, 



(36) 



where the equality holds if pi = P2. If a ligand is bound 
in phase C of a forward cycle (with probability pi), a 
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molecule A in the reservoir will have been transformed 
into a molecule B when the cycle is completed. This 
means that the average decrease in the free energy of the 
reservoir per cycle is given by piAfi, which is identical to 
the right-hand side of eq. (36). Consequently, it turns out 
that inequality (36) is consistent with the maximum work 
theorem of thermodynamics, and the maximum work is 
obtained if pi = P2- 

3.2 Backward operation 

The angle 6 of the rotor is decreased extremely slowly 
by the external agent in a pseudo-static, backward op- 
eration of the engine. A cycle starts with phase B in 
which 9 is decreased from 27r to 02, and then phase C 
{02 > 9 > 0i) and phase A {0i > > 0) follow. The 
work delivered by the engine in each phase is obtained 
from the corresponding work in a forward operation by 
changing its sign. Consequently, the work in a cycle is 
given by 



if a ligand is not bound in phase C, and by 



P P2 



(37) 



(38) 



if a ligand is bound in phase C. 

Since a ligand is bound with probability P2 in phase C 
of a backward operation, the average work delivered by 
the engine per cycle turns out to be 

Wbw = (1 -p2)W^bw +P2WbV = -P2AH - ^f{p2,Pl), 

(39) 

where function /(.t, y) has been defined in eq. (35). Since 
f{x,y) > as remarked earlier, it is concluded from 
eq. (39) that 



(40) 



where the equality holds if pi = p2- 

When the engine is operated backwards, it works as a 
generator in the sense that fuel molecules A are gener- 
ated from waste molecules B. As a result, the free energy 
of the reservoir is increased. The right-hand side of in- 
equality (40) represents the average increase in the free 
energy of the reservoir per cycle, while the left-hand side 
is the work provided by the external agent. Therefore, 
this inequality tells that the minimum work necessary 
for converting molecules B into A is equal to the change 
in the free energy of the reservoir, which is consistent 
with the minimum work theorem of thermodynamics. 

3.3 Reversibility 

It is evident that a single cycle of a forward or back- 
ward operation discussed above cannot be reversible un- 
less pi = P2 =^ 1. Suppose, for example, that a ligand 
is bound in phase C of a forward cycle and the engine 
is operated backwards after the cycle is completed. The 
engine may or may not carry a ligand in phase C of the 
backward cycle, indicating that individual cycles are not 
reversible. 



A cycle may be reversible, however, in the following 
statistical sense. Considering the average over many cy- 
cles, the probability that the engine carries a ligand in 
phase C is pi in a forward cycle and P2 in a backward 
cycle. Therefore, if the condition pi = p2 is satisfied, the 
statistical properties of the whole system (including the 
engine;, the reservoir of ligands, the heat bath and the 
external agent) at a given value of ^ in a forward cycle 
is exactly the same as that in a backward cycle. This 
condition of reversibility coincides with the condition for 
equalities in eqs. (36) and (40) to hold. This observation 
is consistent with thermodynamics: the maximum work 
obtained or the minimum work needed to provide in a 
process is achieved when the process is reversible. 

It should be remarked that, even in the statistical 
sense, pseudo-static operations of the engine discussed 
above are not reversible if pi ^ P2- The irreversibility 
arises when the phase of operation changes from C to B 
in a forward process, for example: The probability that a 
ligand is bound to the engine is pi iov = 62 — ^ (phase 
C), and right after 9 is increased across 02 it must be 
relaxed irreversibly towards p2, the probability in equi- 
librium iov = 02 + Q (phase i?). In a backward cycle 
the irreversibility comes out upon passing from phase C 
to phase A. This kind of irreversibility cannot be elimi- 
nated no matter how slowly is varied. The same kind 
of irreversibility is discussed in ref. 17, and a similar ir- 
reversibility is pointed out^* to result from the contact 
of a small system with a heat bath. 

4. Mechanism of Energy Transduction 

The present model of a molecular chemical engine is 
simple enough to find out how the engine converts the 
free energy (2) released by the chemical reaction into 
mechanical work. The torque exerted by the engine has 
contributions from the conformational free energy Vo(0) 
of the engine without a ligand and from the effective 
interaction e{0) between the engine and a bound ligand. 
It is evident that the e is responsible for net mechanical 
work in one cycle of operation, because the contribution 
from Vq cancels out over a cycle due to its periodicity (8). 
However, the interaction €(0) appears to have nothing to 
do with the chemical potentials 11 a and /i^ of reactant 
and product molecules, which define the decrease A/z 
in the free energy of the reservoir. Then, how the work 
delivered by the engine may depend on /i/i and /ig? We 
will see that the chemical potentials affect the "timings" 
of switching "on" and "off" the interaction e{0), and that 
the ammmt of work delivered by the engine during a cycle 
depends on these timings. 

We first consider an idealized situation in which ha 
lies well above ei ~ e{0i) and jis well above £2 = ^(^2) 
(see Fig. 3), and k-^T is much smaller than jiA — ei and 
Ms ~ £2- In this circumstance we have pi 1 and p2 — 1, 
and observe that pa{0) increases from zero to one in a 
narrow region around 0a defined by e(^^) = as is 
increased from to ^1, according to eq. (15). Similarly, 
Pb{0) decreases from one to zero around angle 0b defined 
by e{0B) = Mb as is increased from 02 to 2tt. This means 
that the engine catches a molecule A when 6 reaches 
around 9a in a forward cycle, and keeps carrying the 



J. Phys. Soc. Jpn. 

ligand afterwards until 9 reaches 9b, where the hgand is 
released as a molecule B. Hence, the interaction e(^) is 
switched on around 6 = 6a and switched off around 6 = 
9b in this situation. Consequently, the work delivered by 
the engine is estimated to be 

w^^lA-m, (41) 

because the torque due to e{9) is given by —Ae/d9 [the 
second term in eq. (11)]. This result is consistent with 
the rigorous one (30) with pi ~ 1 and p2 — 1- In this 
way, we see that the timings of binding and unbinding 
of a ligand molecule determined by ha and hb control 
the work delivered by the engine: if /i/i is increased, for 
example, then Oa decreases and hence the range of 9 over 
which e(^) is on increases; this results in the increase 
in the work delivered by the engine, since e(^) in this 
interval yields positive torque. 

In the above argument, the temperature is assumed 
to be low enough that the switching of the interaction 
e{9) occur sharply around 9 = 9a and 9 = 9b- If the 
temperature is not so low the switching occur gradually: 
for example, binding and unbinding of molecule A are 
repeated many times over a certain range of 9 around 
9a- However, we can say that e(^) is switched on at ^ = 
9a on average^ because probability pA given by eq. (15) 
regarded as a function of e(9) is "antisymmetric" about 
€ = ^A'- transformation e — ha ^ ^J'A — e results in pA — > 
1 ^ PA- Similarly, it can bo said that e{9) is switched off 
at on average as long as p2 — 1- Hence, the argument 
leading to eq. (41) is still valid even in the case of gradual 
switching, if the timings of switching are understood as 
the averaged ones discussed here. 

In the idealized situation considered above, the engine 
works with high efficiency for arbitrary choices of /i^ and 
Hb 8.S far as conditions pi ~ 1 and p2 — 1 are satisfied: 
the decrease in the free energy (2) of the reservoir is 
almost perfectly converted into mechanical work {W ~ 
A/i). If either of these conditions is not satisfied, the 
efficiency is not always good. For example, suppose that 
fj,B lies well below 62 (p2 — 0), and 11 a well above ci (pi ~ 
1). In a forward cycle under this extreme circumstance, 
the engine catch a molecule A at 9 = 9a on average (since 
Pi ~ 1), and releases a molecule B immediately after 
the process enters phase B (since p2 — 0). Therefore, 
potential e(^) is on for < ^ < ^2 on average, and the 
work delivered in a cycle is estimated as 

Wc^HA- 62, (42) 

which is independent of /xb and less than A/x. This re- 
sult is consistent with the rigorous one given by (30) 
with pi 1 and P2 — exp[/3(/iB — ^2)], where the latter 
expression is obtained from eq. (16) under the condi- 
tion that p2 <C 1. According to eq. (4), A^ can become 
large indefinitely by letting the density of molecule B 
in the reservoir become vanishingly small {hb ^ —00)- 
However, the work remains unaltered and the efficiency 
W/Afj, goes to zero in this limit. A similar, qualitative 
argument is given in Ref. 10. 

Another extreme case of interest may be the situation 
in which conditions 

ei > /iA > Ms > £2, (43) 
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Pi <C 1 and p2 — 1 are satisfied. Such a situation can 
be realized by lowering the density of fuel molecule in 
the reservoir. In a forward cycle under this circumstance, 
a fuel molecule A can be bound to the engine only in 
the vicinity of ^ = ^1 (with very small probability). If 
the engine carries a ligand when the cycle proceeds to 
phase C, a forward rcac;tion A ^ B takes place and a 
molecule B will be released around 9 = 9b (since p2 — 1). 
Therefore, the work delivered in this cycle will be 

I^~ei-/is, (44) 

which is larger than A/i = /iyi— /is according to condition 
(43). This result is again consistent with the rigorous one 
given by (30) with pi -C 1 and p2 — and provides an 
example in which work larger than the maximum work 
predicted by thermodynamics can be obtained in a sin- 
gle cycle of forward operation [the second inequality in 
eq. (32)]. 

5. Concluding Remarks 

We have proposed a simple, idealized model for a 
molecular chemical engine, and analyzed pseudo-static 
processes associated with it. The work delivered by the 
engine in a forward operation (§ 3.1) and the work nec- 
essary to operate the engine backwards (§ 3.2) have 
been calculated analytically. The statistical averages 
of these works are found to be consistent with the 
maximum/minimum work theorem of thermodynamics, 
whereas the work for an individual cycle of operation is 
not necessarily so. We have also clarified (§ 4) how the 
densities of fuel and product molecules in the reservoir, 
which determine the free energy released by the chemical 
reaction, can affect the work delivered by the engine. 

In the present work, we have restricted our attention to 
pseudo-static processes controlled by an external agent. 
It is possible to make the model engine constructed here 
run autonomously (i.e.. without the external agent) and 
do work against a load torque applied to the rotor. The 
angular velocity of the rotor, the thermodynamic effi- 
ciency and other quantities characterizing the perfor- 
mance of the engine are under investigation. 

The recent experiment on Fi-ATPase,^^ a biological 
rotary motor, has clearly shown that it consumes three 
ATP molecules per turn of the rotor (the 7 subunit) in 
a forward operation, while a forced, reversal rotation of 
the rotor produces about 2.3 ATP molecules per turn 
on average. Although the experimental situation cannot 
be considered as pseudo-static, it would be tempting to 
compare it with the present analysis. The experimen- 
tal results mentioned above indicate that pi ~ 1 and 
P2 — 0.77, since Fi-ATPase contains three engines called 
the /? subunits and each engine is thought to be responsi- 
ble for 120° rotation of the rotor. In our model, probabil- 
ity p2 depends on the concentration of product molecules 
in the reservoir. In the case of Fi-ATPase, the product 
molecules are adenosine diphosphate (ADP) and inor- 
ganic phosphate (Pi). The physiological concentrations 
of ATP and Pi were chosen in the experiment. If our 
simple model captures the essence of energy transduc- 
tion in Fi-ATPase, it is predicted that different choices 
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of ADP and Pi concentrations will result in clifFereiit val- 
ues of P2 (mechano-chemical coupling efficiency); p2 will 
increase if these concentrations are increased. 

Finally, wc mention a few theoretical works that have 
closely related with the present work. A quasi-static pro- 
cess of a macroscopic chemical engine is discussed in 
ref. 19. In refs. 17,20 models of molecular motors driven 
by two reservers containing the same kind of particles 
with different densities; the decrease in the total free en- 
ergy of the reservoirs resulting from the transfer of par- 
ticles from the high-density reservoir to the low-density 
one through the motor is converted into mechanical work. 
The Carnot cycle for a microscopic heat engine is ana- 
lyzed in ref. 21. 
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